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We calculate the connected insertions of the nucleon three-point function to study the first few 
moments of the unpolarized structure functions of the nucleon. (The disconnected insertions are 
discussed elsewhere in these proceedings). The calculation employs the CP-PACS/JLQCD 2+1 
dynamical clover fermions on a 16^ x 32 lattice with lattice spacing 0=0.1219 fm. The sequential 
source technique, using non-zero and zero momentum point nucleon field as the secondary source, 
is applied enabling a study of different currents at various momentum transfer. 
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1. Introduction 

The calculation of moments of the light cone parton distributions gives detailed insight about 
the structure of the nucleon. The spin average (unpolarized) structure functions Fi and F2 give 
information about the density of quarks and gluons in the nucleon, which leads to an understanding 
of its structure. Structure functions of hadrons can be calculated only non-perturbatively. Using the 
operator product expansion, it is possible to calculate moments of the structure function. Moments 
of the quark density distribution are related to the following matrix elements of twist-2 operators 



(x"-i) = i(P5|(^)"-V7{M.'o^2--'^M4V^I^^) (1-1) 

Here { ... } indicates symmetrization and D = ^{D — D). The most promising method known 
to calculate and understand the quark and gluon structure (i.e. to calculate these matrix elements) 
is lattice QCD. Lots of calculations have been done in this area in the quenched approximation 
% |, |, [T|, P, III, [ll], where the vacuum polarization effects are neglected. It is important to 



calculate various physical quantities for full QCD [|3|, |, ^ |ri|] (including vacuum polarization 
effects) and compare these with the quenched ones. 

We present a calculation for full QCD with relatively heavy pion mass (~600 MeV). In this 
calculation we will neglect gluonic operators. We are planning to repeat the same analysis for 
"almost physical" pion mass (~190 MeV). 

2. The Lattice approach 

Matrix elements given in Eq. ( |1.1[ ) are calculated by evaluating the connected and discon- 
nected insertions as shown in Fig.l. To calculate the nucleon matrix elements we first compute 
two- and three-point correlation functions defined by 

Gr{t2,p) = T.^a,p{Np{t2,p)Na{to,p)), (2.1) 

Grit2,h,P,^) = Y.^a,p{Npit2,p)ff{^,...^,„^,}{h)Na{to,p)) (2.2) 
a.p 

where is the proton interpolating field, given by 



Nait2,p) = L e-'PhatcV^'^\x){Y^"^\x)Crsv'-'^\x)) (2.3) 

x,a,h,c 

C is the charge conjugation matrix, F is the polarization matrix, and 

%....A:„^,} = W{a:,'^M2-'^M„+i}'/^- (2-4) 

To calculate moments of the structure function one needs to take ratios (R) of three-point to 
two-point functions at zero momentum transfer. These ratios are proportional to the desired matrix 
elements. 
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In this work, connected insertions are calculated using the sequential source technique [16] 
with fixed sink momentum. This technique fixes the nucleon and its momentum, but leaves the 
spatial momentum transfer free. 

There are some systematic errors we need to take into account. The most prominent systematic 
eiTors here are operator mixing, contribution from the excited states, finite size effects, discretiza- 
tion errors, and contributions from the disconnected insertions. (The disconnected insertions are 
calculated elsewhere 



3. Simulation parameters 

We perform our full QCD calculation for clover fermions with C5vi'=l-7610 and K"5=0. 13760 
on a 16^ X 32 lattice at lattice spacing (3=0.1219 fm. Configurations were taken from the CP- 
PACS/JLQCD collaboration [^, |l^]. For chiral extrapolation, we run at three different hopping pa- 
rameters, K = 0.13760 ~800 MeV), 0.13800 (m^^ ~700 MeV), and 0.13825 {ntj, ~600 MeV). 
The calculation of nucleon matrix elements, in particular derivative operators, requires high statis- 
tics. For our calculation, we used 800, 810, and 813 configurations respectively with positive and 
negative sink momenta (±1,0,0). We chose polarization along the z-direction (for the calculation 
of magnetic moment and angular momentum). 



4. Results 

Our results for the unpolarized structure functions are listed in Table 1. These are bare (not 
renormalized) quantities. Results are tabulated for a single source position. We are further plan- 
ning to take three different source positions (0,0,0,0), (8,8,8,8) and (16,8,8,16), expecting that this 
decreases statistical errors by a factor of \/3- Calculations for the remaining two masses are in the 
pipeline. 
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Operators: 



moments 


Components 


sink momenta 
P 


K-=0. 13760 


<X>u 


fc>{44} - 3(0{11} +fc>{22} + 0{33}j 
0{44}-5(0{ll}+0{22} + 0{33}) 
0{41} 


/'A A A\ 

(±1,0,0) 
(±1,0,0) 


A o.Cin i 1 (^\ 
U.jy /(Id; 

0.436(32) 

0.395(59) 




0{41 1 } - 5 (0{422} + 0{433} ) 


(±1,0,0) 


0.082(31) 


<x>f 
<x>f 


©{44} -5(0{ll}+0{22} + 0{33}) 
©{44} -i(0{ll}+0{22} + 0{33}) 
0{41} 


(0,0,0) 
(±1,0,0) 
(±1,0,0) 


0.1799(8) 
0.191(14) 
0.165(28) 




®{411} - ^(©{422} +0(433} ) 


(±1,0,0) 


0.036(16) 



Table 1: Table 1 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 



f \ 



\ 



K = 0.13760 



sink time 



(Red: non-zero sinli momentum) 



<x> = 0.436(0.032) <-[8:ll] 



(Blue: zero sink momentum) 
<x> = 0.397(0.016) <-[6:10] 



# configs = 800 



2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 

t 



Figure 2: First moment for up quark: (044 — | (On +022 + 033)) 
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Figure 3: First moment for up quark: (O41) 
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Figure 4: First moment for down quark: {O44 — \{Oii + O22 + O^j,)) 
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Figure 5: First moment for down quark: {O41) 
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Figure 6: Second moment for up quark: (0411 — 5(0422 + C433)) 
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Figure 7: Second moment for down quark: (0411 — 5(0422 + 0433)) 



5. Conclusion 

We have calculated (x) and (x^) for the unpolarized structure function. The results so far are 
encouraging. In order to compare our results to the experimental ones, we have to renormalize 
them, taking into account operator mixing. 

For the larger lattice, we are planning to include smearing at both sink and source to reduce 
contamination from the excited states. We may incorporate operator improvement to reduce the 
discretization effect. The large lattice size will reduce finite size effects significantly. 

We are also working on electric and magnetic form factors for various momentum transfer and 
quark angular momentum. 
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